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خحطر خشکسالی تهدیدکننده و مخل سیستم‌های کشاورزی» محیط‌زیستی, اجتماعی و اقتصادی حوضه 
آبخیز زاینده‌رود است. پایش خشکسالی با استفاده از تصاویر ماهواره‌ای می‌تواند نمایانگر شدت و گستره 
خشکسالی در مناطق با کمبود داده بارش هواشناسی بوده و کاستی مکانی و زمانی آن را جبران نماید. 
کارایی شاخحص‌های پوشش گیاهی ماهواره‌ای مانند شاحص تفاوت پوشش گیاهی نرمال (1۷1) در 
تعیین میزان و شدت خشکسالی روشن است؛ اما ازآنجاکه واکنش پوشش گیاهی به حشکسالی نسبت به 
زمان واقعی تأخیر دارده پایش پاسخ پوشش گیاهی به خشکسالی هواشناسی بر اساس داده‌های سری‌زمانی 
ماهواره‌ای به برنامه‌ریزی کاهش خطرات خشکسالی کمک شایانی می‌کند. در این پژوهش به بررسی 
مکانی-زمانی خشکسالی منطقه‌ای حوضه آبخیز زاینده‌رود طی سال‌های ۲۰۰۳ تا ۲۰۱۶ پرداخته شد. در 
این راستاء با بهره‌گیری از ۸۶ تصویر سری‌زمانی شاخص 12۷ محصولات ماهواره مودیس از فصل 
رشد منطقه» شاخص وضعیت پوشش گیاهی ([/)۷) که یک شاحص نرمال‌شده جهت ارزیابی خشکسالی 
است. محاسبه گردید. سپس پوشش هر طبقه از شدت خشکسالی حاصل از این تصاویر برای دوره زمانی 


پایش‌شده نعیین گردید. شاخحص بارش استاندارد ([۳) مربوط به ۳۹ اتسششکام حهت ارزیابی شاخحص 
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ماهواره‌ای محاسبه شده و ضریب همبستگی این شاخص با میانگین شاخص ۷6 حاصل ارزیابی گردید. 
میزان تأخیر خشکسالی هواشناسی به پوشش گیاهی با بررسی همبستگی شاخص میانگین ۱۷۷ با 
میانگین ٩۳]‏ و نمودارهای تغییرات داده‌های آ/)۷ با داده‌های 8۳1 و 10۷ ارزیابی شد. نتایج پایش 
مکانی-زمانی نشان داد که قسمت‌های مرکز رو به شمال حوضه آبخیز زاینده‌رود. شرق و جنوب شرقی آن 
(محل قرارگیری تالاب گاوخونی) در طول دوره مطالعه به‌طور مداوم با خشکسالی شدید روبه‌رو بوده و 
در سال‌های ۲۰۰۸ و ۲۰۱۶ این خشکسالی بسیار شدید شده است. همچنین نتایج حاکی از تأخیر زمانی 1 
ماهه شاخحص خشکسالی هواشناسی با شاخص خشکسالی ماهواره‌ای ۷ و معنی‌داری ۰/٩‏ شاخص 
1 با شاخص ]۷ است که نشان می‌دهد می‌توان از شاخص‌های ماهواره‌ای در نبود شاخص‌های 
هواشناسی برای بررسی خشکسالی حوضه آبخیز مورد مطالعه استفاده نمود. 
کلیدواژه‌ها: پوشش گیاهی تأخیر پاسخ» 1۷101915-31۷ همبستگی پیرسون. 
۱. مقدمه 
خحشکسالی یکی از مخرب‌ترین مشکلات محیط‌زیستی است که در حوضه آبخیز زاینده‌رود رخ می‌دهد 
(میراحسنی و همکاران. ۱۳۹۵). خشکسالی یک بلای طبیعی نامحسوس است که پوشش گیاهی, زمین‌های 
کشاورزی و منابع طبیعی را در معرض خطر جدی قرار می‌دهد. این پدیده طبیعی و پرتکران مناطق مرطوب و 
خشک را تحت تأثیر قرار می‌دهد (مدرس و سیلوا؛ ۲۰۰۷) و عواقب نامطلوب وسیعی بر شرایط اجتماعی و 
اقتصادی مردم مناطق در معرض خشکسالی می‌گذارد که شامل تأثیر بر دسترسی آب و کیفیت و سلامت اکوسیستم 
است (اسوبودا" و همکاران» ۲۰۱۲). خشکسالی همواره باعث بسیاری از تحولات اقتصادی, جنگ‌ها قحطی و 
مهاجرت‌ها گردیده و به نظر می‌رسد در دهه‌های اخیر فراوانی اين پدیده از نظر شدت. مدت. فضای تحت پوشش. 
تلفات جانی. خسارات اقتصادی و اثرات اجتماعی درازمدت بیشتر از سایر بلایای طبیعی مانند سیل و زمین‌لرزه بوده 
است (خداقلی و همکاران. ۱۳۹۳). درک ویژگی‌های خشکسالی یک گام ضروری برای فراهم آوردن اقدامات 
کارآمد کاهش مشکلات ناشی از خشکسالی بوده و عنصر اساسی در برنامه‌های صحیح مدیریت خشکسالی است 
(مدرس و سیلو؛ ۲۰۰۷). شناخت این پدیده ویژگی‌ها و خصائص آن. آگاهی از امکانات. توانمندی‌ها و 
استعدادهای مناطق ازیک‌طرف و ارائه راهکارهای کوتاه و بلندمدت از طرف دیگر, مدیران منابع آب‌های سطحی را 
بر آن می‌دارد تا به‌منظور کاهش خسارات خشکسالی آمادگی لازم را به وجود آورند. توزیع منطقه‌ای وضعیت 
خشکسالی نسبت به بررسی ایستگاهی آن. سبب درک به‌مراتب کامل‌تری از این پدیده می‌شود؛ بنابراین. آگاهی از 
شدت. تداوم و مساحت گستره‌های درگیر با خشکسالی و احتمال وقوع آن در دوره باز گشست‌های معین می‌تواند در 
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مدیریت منابع آب و در شرایطی که این پدیده منابع حیاتی (آب. خاک گیاه) را در خطر انهدام و نابودی قرار 
داده است» بسیار مفید باشد و مدیران مربوطه را در اتخاذ تصمیمات مناسب و سنجیده در این راستا یاری داده» تا 
خسارات ناشی از سال‌های همراه با خشکسالی را به حداقل برسانند (رضایی و همکاران۱۳۹۳). یکی از راه‌مای 
تعدیل خشکسالی, ارزیابی و پایش آن بر اساس شاخص‌هایی است که بتوان میزان شدت و تداوم آن را در یک 
منطقه تعیین نمود. انواع شاخص‌های خحشکسالی مبتنی بر ایستگاه با استفاده از اندازه‌گیری دما و بارش مانند شاخحص 
بارش استاندارد" )٩۳(‏ توسعه یافته‌اند (مک‌کی " و همکاران. ۱۹۹۳). بر این اساس, شاحص ٩۳1‏ به‌طور موّثر شرایط 
خشکسالی ایستگاه‌های هواشناسی را ارزیابی می‌نماید. شاخص ٩۳1‏ ۱۲ ماهه که در این مطالعه به‌کار رفته است 
نشان‌دهنده الگوهای بارش بلندمدت است که معمولا با جریان‌ها, ذخایر و حتی سطح آب‌های زیرزمینی گره خحورده 
است (2/10, ۲۰۰۷). پایش و ارزیابی پویایی پوشش گیاهی با استفاده از شاحص‌های مشتق از تصاویر ماهواره‌ای 
در شرایط خشکسالی و معمول به کار رفته است. کارایی شاخص‌های خشکسالی ماهواره‌ای مانند شاخحص تفاوت 
پوشش گیاهی نرمال " (0117۷1 (تاکر؛ ۱۹۷۹) در مطالعات ابت شده است (تاکر و چادوری* ۱۹۸۷؛ کوگان" و 
اونگانای. ۱۹۹۸؛ باجگیران و همکاران ۲۰۰۸). بااین حال ازآنجاکه شاحص ۱12۷1 برای پاش خشکسالی شامل دو 
محدودیت اصلی؛ تأخیر زمانی آشکار بین بارندگی و پاسخ 9۷1و نیز کم‌بودن تأثیر بارش در فصل رشد (دوره 
تولید دانه گیاهی) بر 12۷1" می‌باشد. همچنین این مساله که ۷1 نمی‌تواند حطاهای ناشی از انعکاس خاک بر 
گیاه را کاهش دهد (تیان" و همکاران ۲۰۱۵ شاحص وضعیت پوشش گیاهی (۷۲) (کوگان. ۱۹۵۵) که نرمال 
شده‌ی شاخحص ]112۷ است در اين پژوهش به کار برده می‌شود. علاوه بر موارد فوق» شاخحص ۷۲ مبتنی 
بر 70۷1 می‌تواند حطاهای احتمالی ناشی از استفاده از شاخحص ]170۷ را به حداقل برساند. شاحص ۷۲ برای هر 
پیکسل بر اساس حداقل و حداکثر مقدار ۷72۷1 در طول زمان جهت ارزیابی نسبی تغییرات در سیکنال ۱۲۷ 
محاسبه می‌گردد. 

گیتلسون (۱۹۹۷ ثابت کردند که »۷ بر اساس شاخحص ۷۲" نیز با چگالی پوشش گیاهی. زیست توده 
و اندازهگیری بازتاب قابل مقایسه بوده و یک شاخحص خوب تأثیر اقلیم در وضعیت و سلامت پوشش گیاهی 
می‌باشد. گانش (2007) شاخحص ۷۲ را در فصل رشد (ماه‌های مارس تا آگوست) برای ایالت تگزاس تهیه 
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و آنالیز همبستگی آن‌ها را با شاخص‌های خشکسالی هواشناسی ]۳۳5 و ]5۳ ۳ ماهه ۰ ماهه و ٩‏ ماهه 
بررسی نمود که نتایج حاکی از همبستگی بیشتر این شاخص با ]5۳ * ماهه داشت. این عدم همبستگی 
شاحص ]۳۳5 با ۷1 نشان می‌دهد که شاخص ۷6۲ بیشتر شبیه به شاخص‌های خشکسالی که شرایط 
رطوبت را مقدم می‌دانند. می‌باشد. ترا اوره" و همکاران (۲۰۱۶ بیان نمودند که نقشه خشکسالی پوشش 
گیاهی بر اساس شاخحص ۷۲ قابلیت پایش خشکسالی را داراست. دوتا؟ (۲۰۱۵) با استفاده از تصاویر 
بلندمدت 12۷1 ماهواره‌ای خشکسالی کشاورزی را از نظر زمانی و مکانی پایش نمودند. مقایسه مقادیر 
با ٩۳1‏ نشان از همبستگی با ضریب ۰/۷۵ داشت. که کارایی این شاحص سنجش ازدور را برای ارزیابی 
خشکسالی کشاورزی پشتیبانی می‌کند. جیائو* و همکاران (۲۰۱0) ابت کردند که ۷۵۲ می‌تواند شروع 
خشکسالی و تأثیر خشکسالی را تشخیص دهد. به بیان آنان مزیت آشکار ۷7 سهولت محاسبه می‌باشد با 
توجه به اين واقعیت که برعکس شاخص ]33 به مشاهده ایستگاه نیاز ندارد. و به عنوان محصول ماهواره‌ای 
از خشکسالی. اطلاعات نزدیک به زمان واقعی در مقیاس جهانی با قدرت تفکیک نسبتا بالا ارائه می‌دهد. 
روش تلفیقی آنان به عنوان شاخحص خشکسالی بهینه بر اساس سنجش‌ازدور برای خشکسالی طولانی-مدت 
در منطقه گرم کوناس در ایالت متحده آمریکا توصیه شد. زامبرانو" و همکاران (۲۰۱0) با استفاده از 
محصولات مودیس" ۷0۳1301 شاخص ۷۲ را محاسبه نموده و با شاخص ٩۳]‏ مقایسه نمودند. نتیجه 
مطالعه آنان همبستگی بالای داده ها و استرس در پوشش گیاهی ناشی از عدم بارندگی در ماه سپتامبر در ماه 
توامبر را نشان داد. 
پژوهش‌هایی در مورد کاربرد ۲۷ در ۷ برای پایش شرایط خشکسالی در ایران گزارش شده‌است؛ از 
جمله باجگیران و همکاران ۲۰۰۸؛ شمسی‌پور و همکاران ۲۰۱۵؛ اورنگی و همکاران ۲۰۱۱؛ موذن زاده و 
همکاران (۱۳۹۱)؛ رضائی بنفشه و همکاران (۱۳۹۶) و خسروی و همکاران (۱۳۹۶). باجگیران و همکاران 
(۲۰۰۸. با استفاده از نقشه‌های ۷ روزه شاحص 3۲۷۲ ماهواره 0۸۵۸-۵۷۲118۴ به بررسی کارایی 
شاخعص ۷۲ در خشکسالی پوشش گیاهی شمال‌غربی ایران در دوره ۵ ساله و دوره رویشی ۲ ماهه 
پرداختند. همبستگی مقادیر ۷۵۲ با 119*۷1 و مقادیر بارندگی بررسی شده و تأخیر ۳ ماهه شاخص ۷0۲ به 
شاخص ۱2۷۲ حاصل شد. بارندگی هر ایستگاه همبستگی بهتری با مقادیر ۱17۷۲ داشت. نتایج پژوهش 
آنان حاکی از کارایی شاخحص‌های )۷ و ۳۷۲ جهت پایش خشکسالی بوده. که تغییرات بارندگی را با 
6۵۷6۲10 )ظمناهرنا تعحملعظ 9 
6۵ |11 
۵ 2 
0 3 


مصحتصاحصع2 4 
۹۵۵۵۲0۲۵0 عصم 1۳2 ممتانا ۱۷۱۵06۲216۲650 5 


سال ش* بان 5 لی منطقه‌ای حوضه آبخیز زاینده‌رود بر ... ۵ 


پیشآگاهی خشکسالی برای مدیریت ریسک خشکسالی به خوبی نمایش می‌دهد. اورنگی و همکاران 
(۲۰۱۱) در بخشی از مطالعه خود با استفاده از تصاویر ۱۰ روزه ]۲۷ ماهواره ۸۷۳1۳۶ شاخحص ۷۲ را 
برای استان فارس محاسبه نموده و جهت ارزیابی آن» تصاویر حاصل را با نقشه‌های پوشش گیاهی ماهواره 
7 مقایسه نمودند. شاخص ٩۳]‏ جهت تعیین مقادیر کاهش بارندگی در منطقه به کار رفت. نتایج حاکی 
از آن بود که روش به کار رفته پتانسیل پیش‌آگاهی جهت مدیریت ریسک در مناطق نیمهسخشک مانند فارس 
را دارد. شمسی‌پور و همکاران (۲۰۱۵). با استفاده از نقشه های پوشش گیاهی حاصل از تصاویر ۸۵۷۳1۲ 
به پایش خشکسالی طی یک دوره ٩‏ ساله در فصل بهار در دشت کاشان پرداختند. نتایج نشان داد که 
شاخص‌های حاصل از باندهای حرارتی حساسیت بیشتری به شرایط خشکسالی نسبت به باندهای مریی در 
منطقه دارند. شاخحص‌های باندهای انعکاسی مانند ۷ و ]9۷ همبستگی بهتری با پارامترهای هواشناسی 
نسبت به شاخص‌های حاصل از باندهای حرارتی دارند. 
حوضه آبخیز زاینده‌رود به عنوان یکی از مهم‌ترین آبخیزهای کشور. دارای نقش حیاتی در بوم‌شناسی؛ 
معیشت» شرب گردشگری و صنعت کشور است که بررسی و پایش خشکسالی در آن ضرورت مطالعه را 
روشن می‌نماید. یکی از دلایل مهم استفاده از شاخعص ۷۲ در بررسی خشکسالی در حوضه آبخیز 
زاینده‌رود. بر اساس شرایط ویژه پوشش گیاهی این حوزه بوده است؛ زیرا گستره پراکندگی مکانی و تراکم 
و نوع پوشش گیاهی آن از غرب تا شرق حوضه متفاوت است. به‌طوری‌که بر اساس مطالعه دروگر" و 
همکاران (۲۰۰۱) میزان شاخحص ۱7۷1 از منفی (وسعت زیادی بدون پوشش گیاهی و يا به شکل بوته‌های 
تنک) در شرق حوضه آبخیز تا مقدار ۰/۵ در غرب آن به شکل جنگل زاگرس متغیر می‌باشد و در سال‌های 
اخیر حتی در دوره ترسالی از اين میزان تجاوز نکرده است. همچنین بنا بر مطالعه مشترک سازمان فضایی 
ایران " و دفتر پشتیبانی فضایی منطقه‌ای سازمان ملل " که گزارش آن توسط انصاری و همکاران (۲۰۱۵) ارایه 
شده است. استفاده موثر از اطلاعات فضایی برای پایش بلایای طبیعی و تأثیرات آن در مورد خشکسالی کل 
ایران بررسی گردید. در مطالعه مذکور. شاخص‌های ۲۳۷1 *91 10۲ ۷۲ و "۷۳۲ توسط متخصصین 
کشور پيشنهاد داده شد. و نقشه‌های حاصل از تصاویر مودیس به تفکیک استان تهیه گردید. نقشه‌های 
پراکندگی شدت ۷۲ نشان داد که حوضه آبخیز زاینده رود در استان اصفهان با خشکسالی جزیی» و در 
استان چهارمیحال و بختیاری با خشکسالی متوسط رو به رو است. در حالی که شدت این پراکندگی برای 
۱۵ 
۸2 50۵6۵ صهلصه1۳ 2 
۱۷۵-6۲ 3 
همه ] ددعصاطع1ظ 4 
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شاحص ]7۷ متفاوت بوده و از درجه پایین‌تری برخوردار است. ازآنجاکه ۱۲۷ در بررسی خشکسالی 
دارای محدودیت‌هایی که ذکر گردید می‌باشد. نرمال کردن میزان شاخحص 3۲۷۲ که شاحص پوشش گیاهی 
۵ را تشکیل می‌دهد. در کاستن از خطاهای احتمالی ناشی از میزان پایین 3۳9۷1 در حوضه آبخیز مثر 
است و منجر به کارایی این مدل در منطقه می‌باشد. 

هدف این مطالعه بررسی پراکندگی مکانی خشکسالی پوشش گیاهی و ارزیابی کارایی شاخص پوشش گیاهی 
ماهواره‌ای ۷ با استفاده از پایش سری‌زمانی ماهانه تصاویر *«11 ماهواره مودیس و بررسی تأخیر پوشش 
گیاهی نسبت به خشکسالی هواشناسی با محاسبه رابطه ۷ با 5۳1 و ]۲۷ و همچنین بررسی نمودار 
تغییرات شاخص‌ها در حوضه آبخیز زاینده‌رود طی سال‌های ۲۰۰۳ تا ۲۰۱۶ می‌باشد. در اين مطالعه بررسی 
می‌شود که آیا روابط همبستگی بین شاخحص هواشناسی 5۳1 با شاخص‌های ماهواره‌ای )۷ و ۷۲« 
معنی‌دار است» تا بتوان از شاخحص‌های ماهواره‌ای در نبود شاخص‌های هواشناسی برای مطالعات پایش 
خشکسالی استفاده نمود؟ 


۲ منطقه مورد مطالعه 

منطقه مورد مطالعه حوضه آبخیز زاینده‌رود با در اصل حوضه آبخبز گاوخونی با مساحت 1۱1۸۵,1۵ 
کیلومترمربع در فلات مرکزی ایران با اقلیم خشک و نیمه خشک است. این حوضه بین طول جغرافیایی ۲۵ *۵۰ تا 
٩۳ ۶‏ و عرض جغرافیایی "۱۱ ۳۱ تا ۶۲ ۲۳ قرار گرفته است. رودخانه زاینده‌رود به طول ۳۵۰ کیلومتر در 
امتداد تقریبی غرب-شرق در جریان می‌باشد و از کوه‌های زاگرس در غرب استان اصفهان سرچشمه گرفته و به 
باتلاق گاوخونی در شرق می‌رسد (شکل () این رودخانه آب آبیاری, شرب و صنایع مهم کشور را تامین می‌کند. 
اقلیم بیشتر حوضه دارای اقلیم خشک و نیمه‌خشک می‌باشد. متوسط درجه حرارت ماهانه از ۲ درجه سانتیگراد در 
ژانویه تا ۲۹ درجه سانتیگراد در جولای متغیر است. رشته کوه‌های مرتفع زاگرس در غرب و مناطق پست کویری در 
شرق تغییرات مکانی شدید بارش را از بیش از ۷۰۰ تا ۱۵۰۰ میلیمتر (غالبا به شکل ذوب برف) در بهار تا کمتر از 
۰ میلیمتر به دنبال داشته است. میزان بارش از سالی به سال دیگر تفاوت دارد و خشکی‌های مکرر شناخته شده 
است (خداقلی و همکاران ۱۳۹۳). میانگین بارش ۱۲۰ میلی‌متر می‌باشد. از این میزان بارش در ماه‌های سرد و 
زمستان همراه با حرکت جبهه شرقی از شرق اروپا حادث می‌شود و فقط در بعضی مواقع در اثر بارش کافی رواناب 
مهمی در دشت‌های رسوبی ایجاد می‌گردد (سالمی و راست ۱۳۸۳). تیپ پوشش گیاهی منطقه غالبا شامل 
جنگل‌کاری دست کاشت. درختکاری, تالاب. مخلوط تالاب مرتع خوب. متوسط ضعیف. مخلوط مرتع متوسط 
زمین کشاورزی, مخلوط کشت دیم. مخلوط آیشباغ, مخلوط آیش, جنگل تنک. جنگل متوسط باغات. مخلوط 
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کشت اغ. مخلو ط باغات و پوشش گیاهی مخلو ط می‌باشد (میراحسنی؛ 2 زاینده‌رود بعد از طی مسیر حدود 
۰ کبلومتر به تالاب گاوخونی که طبق کنوانسیون رامسر یکی از تالاب‌های مهم ایران به شمار می‌رود و بخش 
زیادی از آن به دلایل مختلف خحشک شده است. می‌ریزد. تالاب گاوخونی نقش حیاتی در تضمین توسعه پایدار 
ساکنان بخش مرکزی ایران دارد. از سوی دیگر از نقطه نظر حفاظت این تالاب به عنوان یکی از اکوسیستم‌های 
ارزشمند در مرکز ایران محسوب می‌شود. پوشش گیاهی در اطراف تالاب بسیار محدود و عمدتا محدود به گونه‌های 
پوشش گیاهی و آبیاری وجود دارد. علاوه بر اين» در مورد تیه همای ماسه‌ای که در اطراف تالاب واقع شده‌اند 
اکوسیستم تالاب گاوخونی در جلوگیری از گرد و غبار تأمین آب تصفیه آب و حفظ ثبات آنها نقش چشمگیری ایا 
می‌کند. لازم به ذکر است که اخیرا خشکسالی‌های مکرر مشکلات آلودگی و محیط‌زیست را افزایش داده و منجر به 
عواقب اگوار در تالاب گاوخونی شده‌است. همچنین تغییر شرایط اقلیمی این اکوسیستم طبیعی را به یک نمکزار 
تبدیل می‌کند (سلطانی و سرحدی» ۳ این مشکلات باعث افزایش خحطر خحشکسالی در حوضه آبخیز زاینده‌رود 
شده که نشان‌دهنده نیاز مبرم به ارزیابی و پایش خحشکسالی و ایجاد استراتژی‌هایی برای کاهش اثرات خشکسالی 


است. 
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شکل ۱. موقعیت حوضه آبخیز زاینده‌رود در استان اصفهان و پراکنش مکانی ایستگاه‌های بارانسنجی مورد 
مطالعه 
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۳ مواد و روش‌ها 

داده های مورد استفادهبه‌طور خلاصه در اين پژوهش از داده‌های سری زمانی ماهانه بارش ایستگاه‌های باران 
سنجی محدوده حوضه آبخیز زاینده‌رود جهت بررسی خشکسالی هواشناسی و از داده‌های سری زمانی ماهواره‌ای 
تصاویر ماهانه مودیس 31۳۷۲ (۷0113۸3) جهت بررسی منطقه‌ای میزان تأخیر خشکسالی هواشناسی بر 
خشکسالی پوشش گیاهی استفاده شده است. 


داده هواشناسی حهت تعیین شاخص بارش استاندارد ([۳) 

جهت تهیه شاخحص 5۳1 نیز ازآنجاکه حداکثر پراکندگی مکانی ایستگاه‌های هواشناسی (۲۱ ایستگاه) در منطقه و 
حداکثر بازه زمانی دوره آماری بلند-مدت ممکن و موجود با شروع یکسان مدنظر بود از داده‌های سری زمانی ماهانه 
سال‌های ۲۰۰۳ تا ۲۰۱۶ (۱2۶ داده) استفاده شد (رجوع شود به شکل ۱). شاخحص 8۳ برای اندازه‌گیری کمبود 
بارندگی استفاده می‌شود. رطوبت خاک در زمان کوتاه‌تری شرایط غیرنرمال بارندگی را منعکس می‌کند. درحالی که 
آب زیرزمینی» جریان‌های سطحی و ذخاثر آب سطحی در زمان طولانی‌تری نسبت به آن عکس‌العمل نشان می‌دهند. 
ازنظر ریاضی» ٩۳1‏ بر مبنای احتمالات تجمعی بارندگی در یک ایستگاه استوار می‌باشد. تام (۱۹۵۸) دريافت که 


1 
تک عسسد < ( اراق 
رابطه (۱) (0) 2۳1 


در اين رابطه ۶ مشخصه شکل و 0 < ۵1 . و 7 میزان بارش و 0 < 3 , 


تابع گاما است که از رابطه زیر بدست می‌آید: 


رابطه (۲) 

این توزیع چوله به راست با حد پایین صفر مشابه توزیع فراوانی بارش می‌باشد که به‌تدریج فراوانی بارش‌ها 
فزایش و در نقطه‌ای به اوج می‌رسد و سپس کاهش و در نقطه‌ای به صفر می‌رسد برای محاسبه 5۳ لازم است. تابع 
چگالی احتمال توزیع گاما را بر توزیم فراوانی مجموع بارش‌ها (برای هر مقیاس زمانی) در هر ایستگاه برازش دهیم. 
مقادیر ۰8 0 تابع چگالی احتمال توزیع گاما برای هر ایستگاه برای هر مقیاس زمانی مدنظر برآورد می‌شود (مک 
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کی ؛ ۱۹۹۳). ازآنجاکه تابع گاما برای مقادیر صفر تعریف شده نمی‌باشد و توزیع بارندگی ممکن است دارای مقادیر 
صفر باشد. تابع احتمال تجمعی که در بر گیرنده مقادیر صفر هم باشد به‌این صورت می‌شود: 
رابطه (۳) 6( -0+ 4 < ()1] 

دراین رابطه0 احتمال وجود صفر در مقادیر بارش است. اگر 19 شمار صفرها در دوره زمانی بارش باشد. تام" 
(۱۹77) ثابت می‌کند که و از تقسیم10 بر می‌تواند محاسبه شود. تام جدول‌های تابع گامای ناقص را برای تعیین 
احتمال تجمعی() [) به‌کار برد. پس از محاسبه تابع تجمعی کل (26) 3 تغییر شکل هم احتمالی تابع تجمعی گاما 
به متغیر تصادفی نرمال استاندارد 2 با ميانگین صفر و واریانس یک صورت می‌گیرد که این متغیر همان مقادیر ٩۳]‏ 
بوده و یک تبدیل با احتمال یکسان است. چنانچه شمار ایستگاه‌ها زیاد باشد و بخواهيم برای همه مقیاس‌های زمانی 
و برای هر ماه سال چنین نموداری را رسم کنیم. بسیار سخت بوده و به آسانی امکان‌پذیر نیست. بدین منظور 2 یا 
۳1 به آسانی از تقریب ابراموویتز و استگون " بدست می‌آید که احتمالات تجمعی را به متغیرتصادفی نرمال استاندارد 


2 تبدیل می‌کند و به صورت روابط زیر است: 


۱ رت + + وت و رو 
3 2 
رابطه (ع) اب + رل +011 +1 
پرای 0.5 > (6) ۲۲ >0 
۱ کات + و یت زو 2 
3 2 
رابطه (۵) + 1 +0 +1 
راطه () ( 0 
پرای 0.5 > (6) ۲۲ >0 
9 
رابطه ۵ (و) 7 -1.0) 


برای 1.0 > (6) 1 > 0.5 
و ثابت‌های معادله عبارتند از: 


,20.010326 ,6(20.802853 ,15517 6022.5 
8و1 ,0.189269حو1 ,۰432788 1,1 
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از نظرمفهومی 8۳1 نمره سیگمایی یا میزان انحراف استاندارد بالاتر یا پایین‌تر از میانگین را نشان می‌دهد. شکل 
(۲) نشان می‌دهد که در پایه زمانی که مشخصه‌های گاما برآورد می‌شوند. 8۳1 دارای یک توزیع نرمال استاندارد با 
میانگین صفر و واریانس یک خواهد بود (ابراموویتز و استگون ۱۹71۶؛ رضایی و همکاران (۱۹۹۳). با بهره‌گیری از 
روش تشریح شده در بالا به کمک برنامه راینه‌ای 6_اه_ذمه برای تنظیم شد که با آن مقادیر 8۳1 در پایه‌های زمانی 
مختلف برای همه ایستگاه‌ها محاسبه شد. با توجه به اينکه این شاخحص در پایه‌های زمانی مختلف و منتهی به همه 
ماه‌ها می‌تواند محاسبه و استخراجح شود. ولی به دلیل اینکه در این پژوهش تحلیل فراوانی خشکسالی سالیانه مد نظر 
بود. مناسب‌ترین دوره زمانی برای این‌کار به جهت مطابقت با داده‌های شاحص‌های ماهواره‌ای. شاخحص ٩۳‏ در پایه 
زمانی ۱۲ ماهه منتهی به ماه دسامبر بود. بر این اساس دوره زمانی شاحص ]53 در مقیاس زمانی ۱۲ ماهه منتهی به 
دسامبر محاسبه شد و برای تحلیل فراوانی مورد بهره‌گیری قرار گرفت. شاخحص بارش استاندارد ۱۲ ماهمه که در ایین 
مطالعه به‌کار رفته است نشان‌دهنده الگوهای بارش بلندمدت می‌باشد. که معمولاً با جریان‌ها؛ ذخایر و حتی سطح 
آب‌های زیرزمینی گره خورده است (اسووبودا؛ ۲۰۱۲). ازآنجاکه ۳1 استاندارد شده‌است. برای اقلیم‌های خشک و 
مرطوب قابل استفاده و نتایج آن برای دوره موردنظر قابلمقایسه است. شاخص 8۳ بر اساس داده‌های ایستگاهی 
ننایج مطالعه شدت مدت شاحص بارش استاندارد 3۳ بازه ۱۲ ماهه جهت این مطالعه به کار گرفته شد (میراحسنی 
و همکاران۱۳۹۵). 


۴/5۳۴۱ < ۱۱ ۰ 7 


۳-۱ ۰ 5۳۱ > ۱۱ ۰ 5 


۹۳ 
شکل ۲. مقادیر ٩۳1‏ دارای توزیع نرمال استاندارد با میانگین ۰ و واریانس ۱ (ابراموویتز و استگون. ۱۹۶) 


داده‌های ماهواره‌ای جهت پایش پوشش گیاهی 

در مطالعه حاضر تصاویر فصل رویشی سنجنده ماهواره مودیس در حوضه آبخیز زاینده‌رود که از فوربه سال 
۰ عرضه شده است. استفاده گردید. به‌منظور تعیین فصل رشد از پردازش تشخیص تغییرات پوشش گیاهی در 
منطقه موردمطالعه تصاویر ماهانه ۱۷1 استفاده شده است. بر این اساس دوره رویشی در این حوضه از بایان ماه 
آوریل تا ابتدای دسامبر (يا ابتدای مه تا آخر نوامبر) تعیین شد (میراحسنی؛ ۱۳۹۲). 


و56 24 ۸۵۲۵۲۵۷1۲2 1 
2 2 


سال ش* پایش خشکسالی منطقه‌ای حوضه آبخیز زاینده‌رود بر ... ۱۳۱ 


تغییرات شاخص ۷۲01016-01۷۲ ماهواره تر!۱ 
شاعص 1۷1 رتاک ۱۹۷۹) معروف‌ترین شاعص پوشش گیاهی ماهواره‌ای است با توجه به رابطه ۸ به 
دست می‌آید: 
رابطه (۸) رمبون ‏ رم0 ۱۳۷ 


که در معادله مذکور 11 باند مادون قرمز نزدیک و > باند فرمز می‌باشد. که شماره باند در تصاویر ماهواره‌های 
مختلف متفاوت می‌باشد. تصاویر سری زمانی 10115-177۷1 ماهانه (۷10113۸3) مورد استفاده در پژوهش از 
ماهواره ترا و از وب‌سایت شرکت ناسا" بین سالهای ۲۰۰۳ تا ۲۰۱۶ اخذ گردید (۱۶۶ تصویر). شبکه سین " این 
تصاویر ۱ کیلومتر است. تصاویر ميانگین ماهانه ۷17۷1 محصولات ماهواره مودیس جهت محاسبه ۷ به کار رفت. 

شاخص وضعیت پوشش گیاهی (۷) 

شاخص ۷ (کوگان 1۹۹۵) با استفاده از محدوده ۰ تا ۱۰۰ از شاخص ]۱۲2۷ بلندمدت محاسبه و نرمال 
می‌شود. مقادیر پایین نشان‌دهنده تنش پوشش گیاهی و مقدار متوسط نشان‌دهنده شرایط متوسط و مقادیر بالا 
نشان‌دهنده شرایط بهینه و طبیعی می‌باشد؛ ۷0 با رابطه )٩(‏ تعریف می‌شود: 

رابطه )٩(‏ (۱0۷) - مس ۱0۲ / هم ۵۷ - ۱0۷۲ ۷۵۲ 

در رابطه بالا»سس1 ۱۱۳۷ و منم ۱۷2۷ به ترتیب میانگین حداکثر و حداقل مدت طولانی 112۷1 برای هر پیکسل 

است که از تصاویر ماهواره مودیس با استفاده از ترم افزار 51 ۱۱۸ استخراج شده انْست (حدول ۳6 


حدول ۱. طبقه‌بندی خشکسالی »۷ (کوگان. ۱۹۹۵) 
کلاس | بسیار شدید | شدید | نسبتا شدید | متوسط | نرمال 
۵ <10 | <20[ <30 | <40 | >40 


رابطه پارامترهای اقلیمی و شاخص‌های ماهواره‌ای 

به‌منظور درک وضعیت پوشش گیاهی با توجه به خشکسالی رابطه شاحص ٩۳‏ با مقادیر میانگین ماهانه] ۱۷۲۷ 
و میانگین ماهانه )۷ با استفاده ضریب همبستگی پیرسون " که یک متغیر فاصله‌ای است. بررسی می‌گردد (رابطه 
۰ 


1 2 

2 ۸۵ 

3 1۳ 0 

4 ۳6۵۲80۴ 020 
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رابطه (۱۰) (۲ ۵2۲۱ )00۷ < 7 

در این رابطه (۷ ,26) 0۷ کوواریانس دو متغیره است. و 026 0۷ انحراف معیار هر دو متغیر است. همبستگی 
وجود یا عدم وجود رابطه بین دو متغیر را نشان داده در صورت وجود شدت آن را تعیین می‌نماید؛ به عبارت دیگب 
برای تشخیص حضور و اهمیت ضریب همبستگی جداول با امکانات مختلف استفاده می‌شود. در واقم. رابطه بین دو 
متغیر می‌تواند مثبت يا منفی باشد. ضریب همبستگی پیرسون () بین ۱ و 1-متفاوت است. مقادیر نزدیک ۱+ یا ۱- 
تشان م ده که وایستکی این و تفر بیکر اشت: 

در این راستاه روابط پوشش گیاهی ۷ با شاحص خشکسالی ٩۳1‏ از طریق ضریب همبستگی پیرسون تعیین 
می‌شود. سپس مساحت پوشیده شده توسط هر طبقه خشکسالی پوشش گیاهی سالانه بر اساس دوره رویشی حوضه 
آبخیز زاینده‌رود تعیین می‌گردد. 
6 بحث و نتایج 

۶ ۱. شاخص بارش استاندارد ([5) 

نمودار نوسانات سری زمانی بازه ۱۲ ماهه شاحص ٩۳]‏ برای ۲7 ایستگاه حوضه آبخیز زاینده رود طی سال‌های 
۳ تا ۲۰۱۶ در شکل ۳ نشان داده شده است (میراحسنی و همکاران ۱۳۹۲). این نمودار حاکی از آن است سال 
۸ ۰۲۰۰4 ۲۰۱۰ و ۲۰۱۱ با خشکسالی بسیار شدید. و در سال‌های ۲۰۰۵ ۲۰۱۲ و ۲۰۱۳ با خشکسالی نسبتا 


شدید هواشناسی رو به رو بوده است. 


شکل ۳. نمودار نوسانات سری زمانی شاخص ٩۳1‏ حوضه آبخیز زاینده رود طی سال‌های ۲۰۰۳ تا ۲۰۱۶ 
: ۲. شاخص وضعیت پوشش گیاهی ([۷۵ 
پایش مکانی شاخص )۲۷ 


نقشه‌های شاخحص ۷ بر اساس مقادیر میانگین سالانه فصل رشد تهیه گردید. جهت بررسی طبقات فراوانی 


تیاه پایش خحشک؟ لی منطقه‌ای حوضه آبخیز زاینده‌رود بر ... ۳ 


خشکسالی نسبتا شدید. خشکسالی متوسط و نرمال تقسیم شدند. شکل نقشه‌های سری‌زمانی میانگین ۷۷ سللانه 
۶ ماهه بین سال‌های ۲۰۰۳ تا ۲۰۱۶ در حوضه آبخیز زاینده‌رود را نشان می‌دهد. بررسی پراکندگی مکانی 
خشکسالی بر اساس نقشه‌ها نشان می‌دهد که بیشتر حوضه آبخیز در طی دوره مطالعه با خشکسالی شدید و بسیار 
شدیل مواجه بوده است و و قسمت جئوب شرقی دایما با خشکسالی بسیار شدید در گیر بوده است., بیشترین گستره 
خشکسالی بسیار شدید متعلق به سال‌های ۲۰۰۸ و ۲۰۱۶ بوده است به‌طوری که قسمت اعظم حوضه آبخیز را 
فراگرفته بوده است. نقشه‌ها همچنین نشان می‌دهد که شروع خنده خشکسالی بسیار شدید هواشناسی در سال ۲۰۰۷ 
(نمودار شکل ۳ قسمت ۱-۳ ) هنوز در پوشش گیاهی حوضه آبخیز نمود نکرده است. همچنین نقشه ها بیانگر آنند 
که حوضه آبخیز بهترین حالت مصون‌بودن از خشکسالی بسیار شدید را فقط در سال ۲۰۱۳ تجربه کرده که پس از 
آن در سال ۲۰۱۶ دوباره شدت گرفته است. متاسفانه طی دوره مطالعه تنها بخش‌های بسیار کوچکی از مرکز حوضه 


آبخیز دچار خشکسالی نبوده است که این وجوب توجه جدی به این مخاطره مهم را نشان می‌دهد. 


شماره بیست و چهارم 


۳3 ۳3 101 13 


راهنمای نقشه: خشکسالی خبلی شدید. خشکسالی شدید. خشکسالی نسبتا شدید. خشکسالی متوسط و نرمال 


۰ سس 
شکل ؛. نقشه میانگین »۷ سالانه فصل رشد در حوضه آبخیز زاینده‌رود (۲۰۰۳-۲۰۱۶) 


پایش زمانی شاخص )۲۷ 

شاخحص ۷۲۲ بر اساس مقادیر میانگین ماهانه ۸۶ ماهه فصل رشد (۷ ماه رشد طی ۱۲ سال» بر اساس رابطه (8) 
و مقادیر جداول ویژگی ۱۷1 با استفاده از نرم افزار 103 ۸6615 و ۳۷15.1 محاسبه و بر اساس جدول ۱ 
(قسمت ۲-۲) طبقه‌بندی شد. نمودار شکل ۵ تغییرات سری زمانی مییانگین شاخحص ۷ ماهانه ۸۶ ماهه در حوضه 
آبخیز زاینده‌رود بین سال‌های ۲۰۰۳ و ۲۰۱۶ را نشان می‌دهد. بر این اساس اغلب موارد خشکسالی در ژانویه فوربه 
مارس که معمولاً فصل معمول بارش در اقالیم خشک و نیمه‌عشک است رخ داده است. خشکسالی بسیار شدید در 
ژانویه ۲۰۰۵ و ۲۰۰۷ و فوربه ۲۰۰۳ تشخیص داده شده است. حشکسالی شدید در ژانوبه ۲۰۰۲ و ۲۰۱۶ و فوربه 
۵ رخ داده است. خشکسالی نسبتا شدید در ژانویه ۲۰۰6 و فوریه ۲۰۰6 ۲۰۰5 ۲۰۰۷ و ۲۰۱۶ دیده می‌شود. 
حشکسالی متوسط در ژانویه 0۲۰۰۳ ۲۰۰۹ ۲۰۱۰ و ۲۰۱۳ فوریه ۲۰۰۸ ۲۰۰۹ ۲۰۱۰ و ۲۰۱۲ و دسامبر ۲۰۰۹ 


۶۹ و ۱۳۰۱ رخ داده است. توضیحات بیشتر در شکل ۸ آورده شده است. 


تیاه پایش خحشک؟ لی منطقه‌ای حوضه آبخیز زاینده‌رود بر ... ۱۵ 


200 200 200 200 200 201 201 201 201 
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شکل ۵. نمودار تغیبرات متوسط ۷61 در فصل رشد در حوضه آبخیز زاینده رود (۲۰۰۳ تا ۲۰۱۶) 


پر اساس نقشه‌های تولید شده مساحت پوشیده مقادیر میانگین سالانه ۷ توسط هر طبقه خشکسالی شامل 
طبقات خشکسالی خیلی شدید. حشکسالی شدید. خشکسالی نسبتا شدید» خشکسالی متوسط و نرمال در حوضه 
آبخیز زاینده‌رود از ماه مه تا نوامبر ۲۰۱۶-۲۰۰۳ تهیه گردید که در نمودارهای شکل 1 و ۷ ارائه شده است. با توجه 
به نمودار ميانگین ماهانه ۷ (شکل *) در طول ۸۶ ماه دوره مطالعه» بیشترین مساحت منطقه متعلق به خشکسالی 
بسیار شدید در ماه آگوست و پس از آن ماه‌های ژوئیه و سپتامبر می‌باشد. وسیعترین دسته متعلق به ماه ژوئن می‌باشد. 
با توجه به نمودار 1 که میزان مساحت هر دسته از ميانگین ۷6 فصل رشد طی سال‌های ۲۰۰۳ تا ۲۰۱۶ را نشان 
می‌دهد. بیانگر آن است که بخش اعظم حوضه آبخیز زاینده‌رود در ماه‌های آگوست و پس از آن جولای و سپتامبر با 
خشکسالی خیلی شدید مواجه بوده است که این مساله از خصوصیات اقاليم خشک و نیمه خشک می‌باشد (مدرس 
و سیلول ۲۰۰۷). 


56۴ با ۷ابال 


شکل 1. پراکندگی مساحت هر طبقه ۷ در طول فصل رشد (۲۰۱-۲۰۰۳) (بر حسب هکتار) 


همچنین نمودارها نشان می‌دهتد که حوضه آبخیز در ماه‌های اکتبر و نوامبر با بیشترین میزان حشکسالی شدید 
مواجه بوده‌است که نشانگر حطر گسترده در این حوضه داشته است. دسته نرمال دارای بیشترین فراوانی مکانی در ماه 
ژوئن است. با توجه به نمودار شکل ۷ که میزان پراکندگی مساحت هر طبقه ۷ سالانه طی سال‌های ۲۰۰۳ تا 


۱۹ جغرافیا و مخاطرات محیطی شماره پیست و چهارم 


۶ را نشان می‌دهد. بیانگر آن است که حوضه آبخیز زاینده‌رود در سال ۲۰۰۸ و بعد از آن ۲۰۱۶ و ۲۰۱۰ با 
خحشکسالی خیلی شدید پوشش گیاهی و در سال‌های ۲۰۱۱ و بعد از آن در سال‌های ۲۰۱۲ ۲۰۱۳ با بیشترین میزان 


خشکسالی شدید پوشش گیاهی مواجه بوده انییشت: دسته نرمال دارای بیشترین فراوانی مکانی در سال ۷ است. 


شکل ۷ پراکندگی مساحت هر طبقه میانگین سالانه [»۷ در (۲۰۱-۲۰۰۳) (برحسب هکتار) 


رابطه شاخص ها و پارامترها بر اساس تغییرات در طول دوره مطالعه 

تغییرات ]۷۲ و ٩۳۲‏ 

همبستگی بین میانگین ٩۳1‏ ماهانه و )۷ متوسط ماهانه در طول دوره مطالعه ۱۲ ساله در حوضه آبخیز 
زاینده‌رود برابر ۰/۷۲۶ حاصل شد. ازآنجاکه با افزايش خشکسالی (کاهش میزان آ۳) میزان شاخص ۷۲ کاهش 
می‌یابد» این همبستگی مثبت می‌باشد. شکل ۸ نشان می‌دهد که تغییرات )۷ میانگین سالانه نسبت به آ۳٩‏ نسبتا 
هماهنگ است؛ اما در سال ۲۰۰۸ با کاهش شدید مواجه شده و سپس افزایش یافته است. همچنین, تغییرات ٩۳]‏ 
میانگین در فصل رشد با آل)۷ از ماه مه تا ژوئن به هم نزدیک شده و از ماه ژوئن تا ماه ژویبه به‌طور منطبق کاهش 
یافته است. کاهش ۷ تا آخر فصل رشد ادامه داشته ولی ٩۳1‏ در ماه آگوست افزايش و سپس در ماه سبتامبر دوباره 


کاهش يافته و پس آن تا پایان فصل رشد افزایش می‌یابد. 


۱۰۰ ۲ 
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۱ ۳ 
اکتبر سپتامبر آگوست جولای ژوئن مه اآوریل ۶ ۲ ۱ ۸ و ا تا ا ۲( 


شکل ۸ نمودارهای تغییرات میانگین ۷ و 5۳1 در فصل رشد (راست) و سالانه (جپ) 


سال ش* بانتن رک لی منطقه‌ای حوضه آبخیز زاینده‌رود بر ... ۱۷ 


جهت بررسی میزان تأخیر زمانی بین مقادیر بلندمدت ماهانه ۸۶ ماهه بین ۷ و ۰٩۳‏ ضریب همبستگی ۷0۲ 


نسبت به 5۳ محاسبه گردید که در جدول ۲ آورده شده است. بر این اساس شاخص ٩۳۲‏ و در ششمین فصل رشد 


بیشترین همبستگی را با شاحص ۷۲ داشت و میزان همبستگی بعد از آن کاهش می‌یابد (جدول ۳). 


جدول ۲. تأخیر زمانی ۷ نسبت به 5۳1 


همبستگی | ماه ۱ | ماه ۲ | ماه ۳ | ماه ۶ | ماه ۵ | ماه " | ماه ۷ | ماه ۸ | بهترین همبستگی 


۵و ۱ ۱۰/۲۹ ۰۳ ۱ ۰/۲ ۱ ۰/۲۹ ۱ ۱۰/۳۶ ۱۰/۳۵ ۱۰/۲۹ ۰/۳۲۲ ماه ششم 


تغییرات ]۷ و ]۲۷ 

شکل ٩‏ نشان می‌دهد که نرخ تغییرات سالانه شاخص ۱2۷۲ بیشتر از شاخص ۷ است» ضمن اینکه تغییرات 
هماهنگ هستند (زیرا مقادیر )۷ نرمال شده هستند). همچنین تغییرات در فصل رشد شاخص‌های ۷ و ۱۱۲۷۲ 
در ماه های اولیه فصل رشد نزدیکتر است و در طول فصل رشد فاصله می‌یابد تابه یک روند مشخص برسد و 


تغییرات سالانه هماهنگ می‌باشد. ضمن اینکه تغییرات هر دو شاخص در فصل رشد کاهشی است. 


توامبر اکتبر سپتامبر آگوست جولای ژوئن مه 
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شکل 4. نمودارهای تغییرات ۷1 و ۷1 فصل رشد (راست) و سالانه (چپ) 


تغییرات ۷۲۷ و ٩1"‏ 

شکل ٩‏ نشان می‌دهد که تغییرات میانگین سالانه شاحص‌های ]۱2۷ و ۳ از سال ۲۰۰۳ تا سال ۲۰۰۸ با روند 
مشابه و فاصله یکسان بوده ولی در سال ۲۰۰۸ با کاهش شدید 9۳1 (افزایش خشکسالی هواشناسی) مواجه شده 
است. بعد از آن تا سال ۲۰۰۹ باوجودی‌که خشکسالی کاهش يافته (افزایش 8 روند تغییرات 17۷1 تا سال 
۰ کاهش است و بعد از آن با فاصله کمتر از قبل از خشکسالی وسیع ۲۰۰۸ افزایش داشته است؛ اما مجددا در 


سال ۲۰۱۶ کاهتن شید تشان می‌دهد: من اینکه با شین کاهش در ستال‌های ۲۰۱۱ نا ۲۰۱۲ ووبه وویست: 
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همچنین شکل تغییرات در فصل رشد بیانگر آن است این دو شاخص در ماه مه با هم تقاطع داشته و با کاهش تا 
پایان فصل رشد بیشترین فاصله را از ٩۳1‏ دارند. ضمن اینکه روند تغییرات کاهشی ]۱۲2۷ منظم و افزایش ٩۳]‏ 
(کاهش خحش؟ لی) نامنظم و به شکل سه بار افزایش و دو بار کاهش است. 


توامبر ‏ اکتبر سپتامبر آگوست جولای ژوئن می 
۲ ۴۰ ۱۴۷۰- 


- ۲۱۴۸۰ ۷۰ 
8۶ - 5 
-۲۱۴۹۰ ۴۴۰ 


-۲۱۵۰۰ ۴۶۰ 


-۲۱۵۱۰ 2-۷۰ 4۵ ۷۱ ۰*۰۴ 


"۰333 


-۲۱۵۲۰ -۰ 


- ۲۱۵۳۰ 2-۰ 


شکل 4. نمودارهای تغییرات ۱۷۲ و ٩۳]‏ سالانه (راست) و فصل رشد (جپ) 


۵ نتیجه گیری 

در پژوهش حاضر رویکرد تحلیل خشکسالی در حوضه آبخیز زاینده‌رود با استفاده از شاحص بارش استاندارد 
5۳1 و ۸۶ ماه تصویر ماهواره‌ای سری زمانی پوشش گیاهی ماهانه 1۷10115-17۷1 جهعت پایش شاخص میانگین 
ماهانه آ/)۷ انجام گرفته و نتایج ارایه گردید. میانگین سالانه هر شاخص با همبستگی ۰/۷ و تأخیر " ماهه شاخص 
آ)۷نسبت به ٩۳‏ مواجه بود که با وجود طیف گسترده بارش در حوضه آبخیز زاینده‌رود (به عنوان مثال بارش 
متوسط ۱۲۸۲/۱۶ میلی‌متر در سرچشمه ایستگاه کوهرنگ در غرب تا ۱۶۹/۸ میلی‌متر در ایستگاه اردستان در شرق 
حوضه بین سال های ۲۰۰۳ و ۲۰۱۶ ارتباط معنی‌داری بین شاخص‌های مورد بررسی مقادیر میانگین ۸۶ ماهه 
پارامترهای )۷ و ٩۳‏ یافت شد. با بررسی نمودار تغییرات سری زمانی شانحص وضعیت پوشش گیاهی ۷ با 
شاخحص هواشناسی آ۳ه به دلیل برقرار بودن روابط می‌توان نتیجه گرفت که شاخحص ۷ را می‌توان جهعت پایش 
و تهیه نقشه شرایط خشکسالی در منطقه مورد مطالعه در مقیاس منطقه‌ای بکار برد با نتبجه مطالعه شمسی‌پور و 
همکاران (۲۰۱۱) در دشت کاشان منطبق می‌باشد» در مطالعه آنان همبستگی شاخص‌های ٩۳]‏ با ۷۲۵۲ برابر ۰/۶۳ 
حاصل شد که بالاتر از همبستگی آن با شاخص‌های ]۷۲1 و 107 بود. همچنین در مطالعه‌های گانش (۲۰۰۷)؛ 
باجگیران و همکاران (۲۰۰۸)؛ دوتا و همکاران (۲۰۱۱)؛ زامبرانو و همکاران (۲۰۱7) نیز همبستگی مناسب بین ۷ 
و ٩۳‏ یافت شد و توانایی استفاده از شاخص ماهواره‌ای را به جای شاخحص خشکسالی هواشناسی در تعیین 


خشکسالی پوشش گیاهی حوضه آبخیز زاینده رود که سوال این پژوهش بود را تایید و آن را تقویت نمود. درنتیجه 
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داده‌های ماهواره‌ای ۷1018 با قدرت تفکیک زمانی بالا برای کشف تغییرات سیمای سرزمین در طول زمان مانند 
پایش سلامت پوشش گیاهی با تحلیل سری‌های زمانی شاخص پوشش گیاهی مفید تشخیص داده شد که در راستای 
نتایج مطالعه انصاری و همکاران (۲۰۱۵)» در کارایی مدل ۷ در حوضه آبخیز زاینده‌رود می‌باشد. در مطالعه 
باعقیده و همکاران (۱۳۹۰)» که تنها سال ۲۰۰۳ با مطالعه حاضر همزمانی داشت. هر دو شاخحص ۱۲۷۲ و ۳1 در 
ماه‌های تابستان کاهش داشتند و شرایط حاکمیت خشکی و کاهش پوشش گیاهی را نشان دادند که با نتایج ایین 
پژوهش (نمودار شکل ۱۰) همخوانی دارد. 

همچنین. نتایج نشان می‌دهد که ویژگی بارش تنها عامل موثر بر مقادیر ۷2۷1 و )۷ نیست. در اکوسیستم‌های 
نیمه حشک. ٩۳‏ اثرات بیشتری بر خشکسالی شدید پوشش گیاهی که در این مطالعه پایش گردید دارند. در اینن 
مطالعه تأخیر " ماهه ٩۳‏ بر پوشش گیاهی را نشان داد. مقایسه تأخیر زمانی پاسخ خشکسالی هواشناسی به پوشش 
گیاهی با استفاده از شاحص‌های ٩۳1‏ و ۷ نشان می‌دهد که حوضه آبخیز زاینده‌رود در سال‌های ۲۰۰۸ ۰۲۰۰۹ 
۰ و ۲۰۱۱ با خشکسالی بسیار شدید و در سال‌های ۲۰۰۳ ۲۰۰۵ ۲۰۱۲ و ۲۰۱۳ با حشکسالی شدید هواشناسی 
رو به رو بوده است. در حالی‌که خحشکسالی بسیار شدید پوشش گیاهی در سال‌های ۲۰۰۳ ۲۰۰۵ و ۲۰۰۷ و 
خشکسالی شید پوشتن گیاهی در سال‌های ۳۰۰۱/۲۰۰۵ و ۲۰۱۶ و خشکسالی تسیا شسلدیل در ۲۰۰۱:۲۶۰۶ 
۷ و ۲۰۱۶ تشخیص داده شده است؛ بنابراین به نظر می‌رسد علاوه بر بارش عوامل دیگری همچون دما تبخیر 
و تعرق. ارتفاع و تامین آب برای آبیاری حوضه همراه با میزان و زمان‌های مسدود کردن و باز کردن دریچه سد 
زاینده‌رود» در سلامت و رشد پوشش گیاهی موثر باشد. همچنین در نظر گرفتن مقادیر برداشت از چاه ها و آبهای 
زیرزمینی؛ با بررسی حشکسالی آب‌شناختی حوضه مورد نیاز است که خارج از محدوده این مطالعه بوده و نیازمند 
تحقیقات بیشتر می‌باشد که برای پژوهش‌های همه جانبه حوضه آبخیز زاینده‌رود در آینده توصیه می‌شود. ازآنجاکه 
حوضه آبخیز زاینده‌رود منابع آب بخش عمده‌ای از فلات خشک مرکزی ایران را تامین می‌کند» بررسی خطر 
خشکسالی اجتماعی و اقتصادی نیز در اين منطقه در مطالعات آینده اهمیت دارد. 

لازم به ذکر است که در اقالیم نیمه خشک بارش بیشتر در زمستان رخ می‌دهد و در تابستان بسیار نادر است؛ و 
به‌طور کلی» روند قابل توجهی در بارش ماهانه بیشتر در طول زمستان و بهار رخ می‌دهد. در حالی‌که بر اساس نتایج 
مطالعه بیشترین مساحت از منطقه که متعلق به دسته خشکسالی بسیار شدید بود در فصل تابستان (به ترتیب ماه‌های 
آگوست؛ ژوئیه و سپتامبر) رخ داده است. ماه‌های مذکور بیشتر فصل رشد در حوضه آبخیز زاینده‌رود را نیز در بر 
می‌گیرد. که محدودکننده و بازدارنده جدی در رشد گیاهان منطقه می‌باشد. از سوییء وابستگی شدید پوشش گیاهی 
تابستانه به بارش بهاره به نوبه خود باعث می‌شود بارش چند ماهه نفوذ بیشتری داشته‌باشد. همچنین از محاسبه 
تأخیر شاخص خشکسالی بر پوشش گیاهی می‌توان نتیجه گرفت که بارش بمثابه عامل غالب در پوبایی خشکسالی 
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پوشش گیاهی در منطقه است. بااین‌حال» این رویکرد باید در اقالیم مختلف آزمایش شود. بازه این مطالعه ماهانه و 
دوره رشد بود تصاویر با قدرت مکانی بالاتر و زمانی کوتاه‌تر (۱۰ روزه) برای افزایش دقت ارزیابی در مطالعات 
بعدی توصیه می‌شود. 

ازآنجاکه نقشه‌های شکل ۶ نشان داد. حوضه آبخیز زاینده‌رود بین سال‌های ۲۰۰۳ تا ۲۰۱۶ در در ماه‌های اکتبر و 
نوامبر با بیشترین میزان خشکسالی شدید مواجه بوده‌است که این مسئله نشانگر حطر گسترده در این حوضه می‌باشد. 
همین‌طور به عنوان نمونه» این نقشه ها نشان داد که تالاب گاوخونی در طی این مدت همواره با خشکسالی خیلی 
شدید روبه رو بوده است و در بین فصول رویشی تنها در ماه نوامبر و در سال‌های ۰۲۰۱۱ ۲۰۱۲ و ۲۰۱۳ از 
خشکسالی خیلی شدید مصون بوده است. لازم به ذکر است. به جز خشکسالی, اين تالاب با چندین مخاطره محبطی 
دیگر از جمله آلودگی. گرد و غبار روبه روست. به‌طوریکه آن را به شکل مخزنی از فلزات سنگین و خطرناک 
تبدیل نموده است. حوضه آبخیز زاینده‌رود طی سالیان انبان و پذیرای پساب کشاورزی گسترده‌ای که در مسیر 
رودخانه صورت می‌گرفته است. و صنایع پرمصرفی از جمله صنایع فولاد که از آب حوضه آبخیز زاینده‌رود استفاده 
می‌نمایند. بوده است. مساله آلودگی به علاوه حشکسالی باعث شده است تا اکوسیستمی که زمانی نه چندان دور 
بسیار غنی پویا و ارزشمند بود و دارای نقش مهم و حیاتی نه تنها برای تفرج محبط‌زیست و اکوسیستم‌های گیاهی و 
جانوری» بلکه برای ساکنین منطقه داشت با نبود آب و برداشت گرد و غبار تهدید جدی برای ساکنان از جمله 
شهرستان ورزنه و حتی استانهای همجوار که در مسیر برداشت تبه های ماسه‌ای و گردو غبار اطراف این تالاب 
هستند باشد. این مشکلات وجوب توجه و برنامه ریزی اضطراری و هرچه سریعتر مسئولان و اجرای سختگیرانه 
قوائین را دو صد چندان روشن می‌نماید. 

در پایان مدل مورد استفاده در پژوهش و استفاده از شاخص‌های ماهواره‌ای در نبود شاخص‌های هواشناسی به 
عنوان راهکار عملی در راستای بهبود کارایی و کاربردی کردن هرچه بیشتر بررسی خشکسالی در اقالیم مشابه اقیلم 


زاینده‌رود و برای کل حوضه‌های آبخیز کشور توصیه می‌شود. 
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